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Abstract 

This paper addresses the problem of reorienting a 
rigid spacecraft from arbitrary initial conditions to 
prescribed final conditions with zero angular veloc- 
ity. The control law analyzed is based on quaternion 
feedback and leaves the user to choose two gains as 
functions of position, angular rate, and time. For 
arbitrary initial states, conditions on the controller 
gains are identified that guarantee global asymptotic 
stability. For the special case of rest-to-rest reori- 
entations, the control law reduces to earlier results 
involving a principal axis rotation. The paper also 
addresses slew rate constraints, both, in terms of the 
two and infinity norms. 

Introduction and Literature Survey 

Spacecraft reorientation problems have been 
treated extensively in the technical literature [1]- 
[8]. Open-loop approaches enable the calculation 
of high-precision solutions that minimize a user- 
prescribed cost index such as fuel consumption or 
maneuver time. However, these approaches usually 
involve iterative procedures and are hence, in most 
cases, computationally too expensive and unreliable 
for on-board implementation. 

Closed-loop or feedback approaches perform only 
noniterative procedures and calculate the current 
control action based on the current state and time. 
Feedback approaches typically perform at best near- 
optimally, and for nonlinear controllers and/or non- 
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linear dynamical systems it is highly nontrivial 
to guarantee that the controller always drives the 
states to the desired position. Additionally, it is 
usually difficult to enforce state constraints, such as 
slew rate limits. 

In practice, slew rate limits often arise from the 
requirement not to exceed the capabilities of a space- 
craft inertial reference unit (IRU) to sense angu- 
lar velocity; slewing too fast can saturate the IRU 
thereby limiting the controller’s ability to track the 
maneuver. Depending on the spacecraft’s archi- 
tecture and the involved on-board instrumentation 
these physical constraints are best captured by two- 
norm or infinity-norm constraints on the angular po- 
sition. 

The general nonlinear feedback control law ana- 
lyzed in this paper has been in service in one form 
or other for many years in numerous applications of 
near-minimum-time satellite reorientation[9]-[ll]. It 
is derived through feedback linearization applied to 
a representation of the dynamical system in terms 
of quaternions in favor of Euler angles. The pa- 
per completely analyzes this controller in terms of 
its physical properties and identifies the conditions 
on the controller gains required to guarantee global 
asymptotic stability and satisfaction of slew rate 
constraints. Thus, the paper provides a frame work 
of conditions within which the control engineer can 
optimize the control gain functions’ remaining de- 
grees of freedom without having to consider global 
stability issues and satisfaction of slew rate limits. 

In [12], Wie and Lu have investigated the con- 
strained slew controller implemented on the XTE 
spacecraft[10]. This controller represents a special 
case of the general controller analyzed in the present 
paper. However, Wie and Lu only investigated rest- 
to-rest maneuvers and assumed perfect execution 
of all control commands without regard to exter- 
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nal perturbations. For this special case, the reorien- 
tation maneuver reduces to a one-degree-of-freedom 
principal axis rotation. The intent of this paper is to 
provide a more comprehensive look at the stability 
of the constrained eigenaxis slew maneuver. 

Problem Formulation 

Consider the attitude dynamics of a rigid space- 
craft. The equations of motion can be represented 
in the following form: 

Juj = — UJ x Juj + u (1) 

1 1 

q = - -w X q + -q A w (2) 

94 = - \^ T q- (3) 

Here, J denotes a constant 3 x 3-inertia matrix, 
ui T = [wi,U2,W3] denotes the angular velocity vec- 
tor, u T = [ui,«2,«3], denotes the control torque 
input vector, and q T = [91,92,93] together with q 4 
represent the angular position through its associated 
quaternion vector [q T ,q^\. 

Euler’s rotational theorem states that the rigid- 
body attitude can be changed from any given orien- 
tation to any other orientation by rotating the body 
about an axis, called the Euler axis, that is fixed, 
both, with respect to the inertial frame as well as 
with respect to the rigid body. For a given reori- 
entation, let e T = [ei,e2,e3] denote a unit vector 
along the Euler axis, and let 9 denote the associ- 
ated rotation angle about the Euler axis. Then the 


quaternion vector [9 T ,94] satisfies 


q 4 = e 4 sin(6/2) 

(4) 

q 2 = e 2 sin(9/ 2) 

(5) 

93 = e 3 sin(9/ 2) 

(6) 

q 4 = cos(9 / 2) 

(7) 


By construction, the quaternion vector [q T , 54] is al- 
ways of norm one, i.e. 

9i + 92 + 93 + 94 = 1- (8) 

The vector [9 T ,94] has no physical meaning unless 
constraint (8) is satisfied. It is important to note 
that Constance of the norm of [q T , 54] is guaranteed 
through the equations of motion (1) — (3), irrespec- 
tive of the chosen control u. Hence, (8) is guaranteed 
to be satisfied (up to the precision of the numerical 
integration) if it is satisfied at a single point, say, at 
the initial time. 

We consider the following problem: 


Definition 1 Given the equations of motion (1)- 
(3) and a positive constant uj max , find a feedback 
control law 

u = f(u,q,q 4 ,t) ( 9 ) 

that drives the states uj, q, q 4 from arbitrary initial 
conditions uj(t 0) E R 3 , q(t 0) E R 3 , 94(^0) £ R with 
q{t 0 ) T q(t 0 ) + rq 4 (t 0 ) 2 = 1 to uj = 0 , q = 0 , q 4 = ± 1 , 
such that the state constraints 

I HO 1 1 2 < u max ( 10 ) 

and/or 

IHOII 00 < U max (11) 

remain satisfied for all times. 

The symbols, 1 1 - 1 1 2 and || • ||oo denote the Euclidean 
and 00 norms in the corresponding finite dimen- 
sional vector space, respectively, i.e. 

IKOIh = \/uii(t ) 2 +Ul 2 (t ) 2 + Ul 3 (t ) 2 

and 

I KOI loo = max{|wi(t)|, MOI, MOI} 


Eigenaxis Rotation Controller 

In the present paper, a controller (9) of the form 
u(uj,q,q 4 ,t) = u> x Juj 

~ k(ui,q,q 4 ,t) Jq 

- c(w, q, q 4 ,t) Juj (12) 

is investigated, where k and c denote scalar func- 
tions of their arguments u>, q, q 4 , and t, respectively. 
This form of the controller is motivated by the de- 
sire to eliminate the cross product term in (1) as 
well as a state feedback. Explicitly, inserting (12) 
back into (1) yields 

uj = —kq — cuj, (13) 

so that (12) can be interpreted also as feedback lin- 
earization if k and c do not explicitly depend on 
states. 

In the following, we will prove a number of useful 
properties of the controller (12). 

Property 1 Consider the dynamical system (1)~ 
(3) with the feedback controller (12), where k and 
c are Lipschitz-bounded with respect to u>, q, q 4 , and 
piecewise continuous with respect to t. Then u>(t), 
q(t), w(0) and g(0) are colhnear if and only if u>( 0) 
and g(0) are colhnear. 
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Proof of Property 1: 

=>■) Assume Lo(t), q(t), w(0) and g(0) are collinear 
for all times. Then, obviously, w(0) and g(0) are 
collinear. 

■$=) Assume w(0) and g(0) are collinear. We first 
show that Lo(t) x q(t) = 0 for all times. We have 

— (lo X q) = — c ui X q — — lo X (lo X q). (14) 

Here we used q x q = 0 and lo x lo = 0. Now con- 
sider the system of ordinary differential equations 
consisting of the equations (2), (3), (13), and (14). 
It can be verified that the right-hand side of this sys- 
tem is Lipschitz-bounded with respect to the states 
q, q 4 , lo, and (lo x q), and independent of time t 
(thus measurable with respect to time t.) Hence, 
whenever a solution to the system (2), (3), (13), 
and (14) exists for given initial conditions, this so- 
lution is determined uniquely (see the Appendix I 
of [13]). As (lo x q) = 0 obviously furnishes a solu- 
tion to equation (14) subject to the initial condition 
(lo x q)( 0) = 0 it is clear that lo x q is guaranteed to 
remain identically zero for all solutions to the sys- 
tem (2), (3), (13), whenever the initial conditions 
are such that w(0) x </(0) = 0. Hence co(t), and q(t) 
are guaranteed to remain collinear whenever w(0) 
and < 2 ( 0 ) are. 

To show that Lo(t) and q(t) are collinear not only 
with respect to each other but also with respect to 
</(0), let us define the new states 

fi(t) = w(f)x{(0), f , 

Q(t) = q(t)xq( 0). 

Then, using the fact that Lo(t) x q(t) is identically 
zero, it can be verified that Q and Q satisfy the 
following differential equations: 

tt(t) = -kQ(t)-ctt(t), 

Q(t) = ±q 4 (t)S2(t). [ J 

Considering the augmented system of differential 
equations consisting of equations (2), (3), (13), and 
(16), and using arguments similar to the ones used 
above to show that Lo ( t ) x q(t) is identically zero, 
we can show that £2 and Q remain identically zero 
whenever they are zero at initial time. This implies 
that Lo(t) and q(t) are collinear to g(0) if w(0) is 
collinear to §(0). 

By virtue of the fact that Lo(t) and q(t) are 
collinear for all times if w(0) and g(0) are, the pre- 
vious result implies that Lo(t) and q(t) are collinear 
also to w(0) if w(0) is collinear to §(0). 


The statement of Property 1 can be strengthened 
by replacing “piecewise continuous with respect to 
time t” through “measurable with respect to time t” . 
Loosely speaking, the set of measurable functions in- 
cludes all piecewise continuous functions, plus some 
functions with an infinite number of discontinuous 
jumps. For a more precise characterization of mea- 
surable functions, see the Appendix I of [13], or Sec- 
tion 3.7 of [14]. Note that by definition, a function 
f(x) is Lipschitz-bounded with respect to x if and 
only if for every x in its domain there is a constant 
L such that \\f(x + Ax) — /(a?) 1 1 2 < L 1 1 Aa? 1 1 2 for 
all sufficiently small Ax. Hence, the condition of 
Lipschitz-boundedness in Property 1 is guaranteed 
to be satisfied if the functions c and k are differen- 
tiable with respect to w, q, q 4 , and have uniformly 
bounded gradients for all fixed t. 

Property 1 shows that w(0), </(()), to(t), and q(t) 
remain collinear for all times if w(0) and g(0) are 
collinear. In this case the body performs a rotation 
about an axis that remains fixed in inertial space. 
The condition that w(0) and g(0) be collinear is triv- 
ially satisfied for a rest-to-rest reorientation where 
w(0) = 0. 

Next we want to investigate what happens if the 
control law in Property 1 is applied to a system 
where w(0) and g(0) are not collinear. Note that 
this situation can arise as a result of external per- 
turbations; it is important that any “out-of-plane” 
component does not grow with time. 

Property 2 Consider ihe system (l)-(3) with the 
control law (12). Now assume there is a constant 
Cmin such that 

c(uj,q,q 4 ,t) > c min y lo, q, q 4 , t. (17) 

Then 

ll w x ( 2 H 2 < ||w(t 0 ) x < 2 (^ 0 ) 1 1 2 • e~ Cm,n (18) 


Proof: We have 

^ (Ik x 9II2) 


(14) (j 1 1 

= 2 (lo x q ) —cuj x q lo x (lo x q) 

T T 

= —2 c (lo x q) (lo x q) — (lo x q) [lo x (lo x </)] 

The second term on the right-hand side is zero as 
(lo x q) T is normal to lo x (lo x q) . Hence 

^ (Ik x g||l) = - 2 c • (w X qf (uo X q) 

= — 2 c | \lo x q\ I 2 
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1 

co< 2 


k 


2 


n -1 


(26) 


With the lower bound (17) on c, this yields 

Ik X l\\l < IK*o) X q(t 0 )\\l ■ e~ 2 Cm ' n 
which implies (18). 


q.e.d. 

Property 2 states that the normal component u> x 
q decays exponentially in norm if and only if the 
constant c msn in (17) satisfies 

Cmin P 0 . ( 19 ) 


£0 < 


£0 < 


do 

T 


2 k 


k 

, K 2 , 


( w tY 


w 0 


C(£2) 


4 (k + C(£2)£2 ) 


do 

2 


(27) 

(28) 


As, by assumption, V approaches d arbitrarily 
closely as time goes to infinity, we can pick a time 
to such that 


d < V < d + €0 it >t o- (29) 


Global Stability Analysis 

Property 3 Consider the equations of motion (1)- 
(3) with the control law (12). Assume that 

A1 ) the gams k and c satisfy the conditions stated 
in Properties 1 and 2, with c m i n > 0, 

A2) k is a positive constant, and 
AS) c is bounded m the sense that for every M £ 
R there is a C(M) £ R such that \c(u>,q,q 4 ,t)\ < 
C(M) as long as ||w ||2 < M. Then, 

PI) lim^oo ||w(2)|| 2 = 0 
P2) lim t ^oo q 4 (t) = +1 or -1, 

PS) lim t ^oo ||g(t)|| 2 = 0. 

Proof: Define the function 

V = 4~k ^ + w 2 + w l) - ?4 (20) 

Differentiating V and inserting (1), (2), (3), (12) 
yields after some calculation 

V = -d- lu t uj < 0. (21) 

z /c 

Assume property P2 is not satished. Since V is 
monotonically decreasing and bounded from below 
this implies that there is a real number d y£ ±1 such 
that V converges to d as time goes to infinity. Let 
do denote the minimum distance of d to ±1, i.e. 

do = min {|rf — 1|, |rf + 1|} > 0. (22) 

Then let us pick positive real numbers £2, wo, and e 0 
such that 

£2 > ^4k(V(t 0 ) + l) 

1 

tiln < — 

- 2 

UJ 0 < \/k do 


The assumption that property P2 is violated (=>• 
q 4 {t) does not go to ±1 =>• ||§(t )||2 does not go to 
zero) implies that we can pick to such that, in addi- 
tion to (29), q(t 0 ) k 0) and we can define 


rp W 0 

1 ~~ 2 (k + C(£2) £2) 

(30) 

T 2 

2 *ll?(*0)ll2 

(31) 


Solving (20) for ui T ui and taking the square root we 
can show immediately that £2 furnishes an upper 
bound on ||w ||2 for all t >t 0 , namely 

IMI 2 = \j4k{V{f) + q 4 (t)) 

< y/4k (y(to) + l) 

< n (32) 

The second line in (32) follows from (29), which im- 
plies V(t 0 ) < v(t), and from \q 4 \ < 1. The third line 
follows from (23). We have to consider two cases, 
namely ||w(t 0 )|| 2 > w 0 and ||w( 2 0 )||2 < w 0 . 

Case 1: Assume 

I k(fo) 1 1 2 > wo- (33) 

From equation (1), with u replaced by the right- 
hand side of (12) we find that for all t >t 0 

Mt)\\ 2 = \\-kq(t)-cco(t)\\ 2 
<k\\q(t)\\ 2 + cMt)\\ 2 
<k + C(£2) £2 (34) 

On the interval from to to to + T\ we hence have 

ui(t) = ui(to) + A ui(t) (35) 

||Aw(t)|| 2 < (* + 0(0)0) Ti 
_ wo 
2 


(23) 

(24) 

(25) 


4 


with 


(36) 



Hence, from (41), (42) 


Hence, for all t £ [to, to + T)] 


ll w (t)|| 2 — IKY + Aw(t)||2 

>|K<o)||2-||Ao;(Oll2 


2 ' 


(37) 


Hence 


V(t 0 +T 1 ) ( = } V(t 0 ) 
< V(t 0 ) 


pto+Ti 


C r r 

— - uj uj at 
2k 


Cmin 

YY 


wo 

Y 


2 

T x 


(28) 

< V (to) — 2 e 0 - (38) 


But, assuming that V(to) satisfies (29), (38) implies 
that V(to + T\) violates (29). Hence case 1 is not 
possible. 

Case 2: Assume 


IHYIh < w 0 . (39) 


\q A - 1| • \q 4 + 1| > (43) 

for all t > to - Noting that 1 54 — 1| • |(/4+l| = \q 4 — 1| = 
q T q, (43) immediately yields 

Il?«ll 2 >y 'it>t 0 . (44) 

From (2) we find now that for all t > t 0 

llslb = II - ^ x q+ 7; 94 ^l | 2 

<f » 

< ^.LO Q < C 1 < C 1 < Clt , o 

<w 0 , (45) 

so that on the interval from t 0 to t 0 + T 2 

q(t)=q(t 0 ) + A q(t), (46) 


If | |oj(i) 1 12 becomes greater or equal wo at some time 
t\ > to, then case 2 reduces to case 1 with to replaced 
by f 1 . Hence, we only have to consider the case where 


|w(t)|| 2 <w 0 Vt>t 0 - 


(40) 


Solving equation (20) for 54 — 1, we find that for all 
t>t 0 


l?4 - 1| 


T 

LU LU 


4 k 

T 

LU LU 


-V-l 


4 k 


> 


(40), (29), (22) 

> 

(25), (27) 

> 


4 k 


- (V - d) - (d + 1) 


— \ V — d| + |d + 1 1 


4 k 


— £0 + do 


do 

Y 


(41) 


and similarly 

l?4 + 1| 


T 

LU LU 


4 k 


4 k 


-V + l 


-(V-d)-(d-l) 


> 

(40), (29), (22) 

> 


T 

LU LU 


4 k 


- \V-d\ + \d-l\ 


-4b-“ +d ° 


(27), (25) d 0 

> — 

2 


(42) 


with 


pto+T2 

1 1 Ag(7) 1 1 2 = || / q(t)dt\\ 2 

Jto 

ft 0+^2 

< / \\ 2 q(t)\\ 2 dt 


< LUo To- 


rn 


Using (46), (47) in equation (1) with u replaced by 
the right-hand side of (12) yields on the interval 
[to, t 0 + T 2 \ 


ui(t) = — k q(t) — c Lu(t) 

= - k q(t 0 ) + A lu (48) 

where 

l|Aw || 2 = \\- k Aq(t) - CLu(t )\\ 2 

< k l|Ag(t)|| 2 + c||w(t)|| 2 

(40), (47) 

< k LUo T 2 + C(LU 0) LU 0 

(32) 

< k LUo T 2 + C(Cl) LUo (49) 

From (48), (49), we see that on the interval t £ 
[to, to + T 2 ] 

Lu(t) = Lu(to) - k q(t 0 ) (t - to) + A Lu(t) (50) 
where 

IIMOH2 

< (k LUo T 2 + C(n) LUo) t 2 
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2 


= w 0 (k Tf + C(n) t 2 ) 

2 

k ||<z(*o)l | 2 

\ 2 


(31) 



w 0 

(44) 


< 

Wo 

(26) 

1 

< 



C(fi) 


k lk(*o)lh 


l. 4a. 

, K 2 , 


C(fi) 


l. 4a 
. 2 , 


(51) 


From (50) we now find 


I Hu, + T 2 ) 1 1 2 

= ||w(H — k q(to) T'j + Aw (to + 72)||2 

> - IHHIh + k | |g(^o ) 1 12 12 

— ||Aw(t 0 + T 2 ) 1 1 2 


(39), (31), (51) 
> 


(24) 

> 


1 


Wo + k Ik(t 0 )|| 2 


2 

k lk(fo )|| 2 



1 

2 


1 

2 


(52) 


But, in light of (24), (52) implies that ||w(to + 
T 2 ) 1 1 2 > wo, which contradicts (40). Hence case 2 
can be excluded, which completes the proof of prop- 
erty P2. 

Property P3 follows immediately from Property 
P2 and q T q + q 4 = 1 (see equation ( 8 )). 

Now assume Property PI is not satisfied. As q 4 — >■ 
±1 (property P2) and V < 0 (see equation (21)) 
this assumption implies that there is a real number 
d > 0 and a time t\ such that ||w(t )||2 > d for all 
times t > t\. With (21) this implies this implies that 
V is unbounded from below. But this contradicts 
Property P2. Hence ||w(t )||2 — >■ 0 as t —>■ 00 . 


q.e.d. 

It is interesting to note that basically the same re- 
sults as those obtained in Property 3 can be derived 
if k is chosen negative instead of positive. Using 

V = —— (uil + uil + uil) + q 4 (53) 

instead of the function V in (20), the proof of this 
proposition is analogous to the proof of Property 3. 
In summary, Property 3 can hence be replaced by 
the following stronger statement: 

Property 4 Consider the equations of motion (1)- 
(3) with the control law (12). Assume that 

Al) the gams k and c satisfy the conditions stated 
in Properties 1 and 2, with c m i n > 0, 

A2) k is a nonzero constant, and 


AS) c is bounded m the sense that for every M £ 
R there is a C(M) £ R such that \c(u>,q,q 4 ,t)\ < 
C(M ) as long as 1 1 w 1 1 2 < M. 

Then, as time t goes to infinity, we have 
PI) | |w(t) 1 12 approaches 0, 

P2) q 4 (t) approaches either +1 or —1, 

PS) |k(t )|| 2 approaches 0. 

Proof: Clear from the remarks above. 


q.e.d. 

Property 4 states that the control law (12) guar- 
antees (under rather mild conditions) convergence 
of the dynamical system (l)-(3) to w = 0, q = 0, 
q 4 = ± — 1. Interestingly, this convergence guaran- 
tee can be made for any constant gain k, both, for 
k > 0 and k < 0. However, Property 4 does not 
provide any help in deciding which sign to pick for 
k to perform a given reorientation most effectively. 
As a rule-of-thumb strategy one may want to apply 
the logic to pick k such that q 4 —>■ +1 if q 4 (to) > 0 , 
and, vice versa. To enable such a choice, the next 
result addresses the question to which value (+1 or 
— 1 ) q 4 can be expected to converge as a function of 
the sign chosen for the gain k. 

Property 5 Assume the gams k and c are constant 
with c < 0 and k y! 0. Then the eigenvalues of 
the system obtained from linearizing the equations 
of motion (l)-(3) with the control law (12) about 
w = 0 , q = 0 , has the following properties: 

PI ) The eigenvalue associated with the q 4 - 
component is always zero. 

P2) If kq 4 > 0, then all eigenvalues associated 
with the w and q-components have negative real 
parts. 

P3) If kq 4 < 0, then three eigenvalues associ- 
ated with the w and q-components have negative real 
parts, and three have positive real parts. 

Proof: Linearizing the equations of motion (l)-(3) 
with the control law ( 12 ) about w = 0, q = 0, leads 
to the linear dynamical system 

6x = A 8x (54) 


where x T = 

kw T , 

Sq T , 

Sq 4 \ 

and 





( ~ c 

0 

0 

-k 

0 

0 

0 

\ 


0 

— c 

0 

0 

—k 

0 

0 



0 

0 

— c 

0 

0 

-k 

0 


A = 

24 

2 

0 

0 

0 

0 

0 

0 



0 

24 

2 

0 

0 

0 

0 

0 



0 

0 

2 ± 

2 

0 

0 

0 

0 



^ 0 

0 

0 

0 

0 

0 

0 

/ 
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Obviously, the eigenvalue A 7 , associated with 5 <74 is 
always zero. Using a determinant relation for par- 
titioned matrices [15], it can be shown that the re- 
maining matrix A, obtained by deleting the last row 
and last column of A, has the eigenvalues 


, , , —2c—\f4c 2 — 8 kq 4 

Ai — A2 — A3 — 

_ _ _ -2 c+ \/4c 2 - 8 k qi 


(56) 

(57) 


Noting that 2 c > \J 4 c 2 — 8 k (74 if kq 4 > 0 , and 
2 c < \/4c 2 -8 k </4 if kq^ < 0 , the statement of 
Property 5 follows immediately from (56), (57). 


q.e.d. 

Clearly, from the results of Property 5 it can be ex- 
pected that under the influence of external pertur- 
bations </4 will always converge to +1 if k is picked 
positive, and, vice versa, 54 will always converge to 
— 1 if A? is picked negative. 

These remarks, however, should be understood 
only as educated guesses. Note that k and c are as- 
sumed constant in Property 5, while c was assumed 
possibly non-constant in Properties 3 and 4. 


Slew Rate Constraint 


To address the slew rate constraint, we will make 
use of the following intuitive logic: The solution x(t) 
of an ordinary differential equation x = f(x,t) is 
guaranteed to satisfy the constraint g(x,t) < 0 if 
the functional dependence of / and g on x and t 
is such that the slope of g goes to zero whenever g 
approaches zero (from below). The associated lem- 
mas summarizing the relevant results are given in 
the Appendix. In the following these results will be 
applied to the slew rate constraints ( 10 ), ( 11 ). 


Two- Norm Constraint 


For the slew rate constraint (10) we have the fol- 
lowing result. 

Property 6 Consider ihe system (l)-(3) with the 
control law (12) and the slew rate constraint (10) 
with 0J m ax > 0. Assume that the functional depen- 
dence of the gams k and c on their arguments lo, q, 
</ 4 , and t is such that 


c > —k 


uj T q 
LO T LO 


(58) 


is satisfied whenever 

ll^l |2 ^ ^ max • ( 59 ) 

Then the slew rate constraint (10), i.e. 

I 1^(0 I I2 C LO max , 

is guaranteed to remain satisfied for all times t >t 0 
if it is satisfied at the initial time to- 

Proof: Noting that 

— (uj T to) = 2 U) T LO 

dt v ' 

= 2 ui T (—kq — clo) 

= — 2 k lo t q — 2 c ui T ui (60) 

the statement of Property 6 follows from Lemma 1 
in the Appendix. 

q.e.d. 

Intuitively, it can be expected that Property 6 re- 
mains satisfied if the inequality sign in (59) is re- 
placed by an equality sign. With the help of 
Lemma (2) it can be shown that this is in fact cor- 
rect if k(u>, q, 54 , t) and c(lo, q, < 74 , t) also satisfy 
the conditions stated in Property 1. We have the 
following result. 

Property 7 Consider the system (l)-(3) with the 
control law (12) and the slew rate constraint (10) 
with 0J m ax > 0. Assume the gams k(u>, q, q 4 , t) and 
c(lo, q, < 74 , t) are Lipschitz-bounded with respect to 
ui, q, </ 4 , and piecewise continuous with respect to t. 
Additionally, assume that k and c are such that 

c>-k^L (61) 

LO 1 LO 

is satisfied whenever 

||^||2 — ^ max • ( 52 ) 

Then the slew rate constraint (10), i.e. 

1 1^(0 1 1 2 ^ ^ max 1 

is guaranteed to remain satisfied for all times t >t 0 
if it is satisfied at the initial time to- 

Proof: Follows from Lemma 2 in the Appendix. 

q.e.d. 

Now assume c and k are such that conditions (58), 
(59) (or (61), (62)) are satisfied. Then it is immedi- 
ately clear that the same conditions remain satisfied, 
if c is divided in half, or if k is multiplied by two. 
In fact, the following general result can be easily 
verified. 
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Remark 1 Assume c and k are such that Prop- 
erty 6 (or Property ? ) is satisfied. Let Ai > 0, 
A 2 > 0 be positive constants. Then, Property 6 (or 
Property ? ) remains satisfied if we replace 

c —>■ c Ai 

k —>■ k A 2 

A 9 

^ max * ^ max ^ 


To see this, note that for ||w|| 2 > oj max 



>1*1 



\^A>-k^L 
IMI1 ~ " T u 


Infinity- Norm Constraint 


From Remark 1 it follows that, by picking c larger 
(or k smaller) than required by the conditions of 
Properties 6 or 7, the maximum attainable slew rate 
is reduced to values smaller than co max . In practice, 
this may represent an undesired performance degra- 
dation, and it may hence be a design goal to sat- 
isfy conditions (58) / (61) only marginally, i.e. with 
equality rather than strict inequality. 

However, simply replacing the inequality signs in 
conditions (58), (61) by equality signs would lead to 
controller gains that maintain ||w(t )||2 at a constant 
value for all times after ||w ||2 = oj max is reached. 
Obviously, this is not desirable. Formally, such an 
approach would also violate condition (17), which is 
required in our proof of global stability (see Prop- 
erty 4). 

A sensible strategy that leads to a controller 
which satisfies Properties 1, 2, 4 of this paper and 
marginally satisfies the slew rate constraint (10), is 
to pick a constant gain k > 0, and to enforce condi- 
tions (58), (61) of Property 6 or 7 with strict equal- 
ity, but to override this constraint (and instead en- 
force (58), (61) with strict inequality) whenever this 
is required by at least one of the conditions of Prop- 
erties 1, 2, 4. It can be easily verified that all condi- 
tions of Properties 1, 2, 4 can be satisfied this way. 
The approach leads to a constant gain k and to a 
piecewise defined, continuous gain function c. 

A convenient simplification of the conditions of 
Properties 6 and 7 is as follows. 

Remark 2 The conditions (58), (59) of Property 6 
(conditions (61), (62)) of Property 1) arc guaranteed 
to be satisfied if the functional dependence of k and 
c on their arguments u>, q, q$, and t is such that 

c= \ k \Ml (63) 

^max 

whenever ||w|| 2 > u max (whenever ||w|| 2 = oj max ). 

For the special case where —kuj T q = \k\ ||w|| 2 ||</|| 2 , 
the right-hand sides of conditions (63), (58) (con- 
ditions (63), (61)) are identical. In this case, 
condition (58) (condition (61)) is satisfied only 
marginally and the slew rate constraint (10) can be 
satisfied with equality. 


In many practical applications, the slew rate con- 
straint has to be enforced in terms of the infinity 
norm (11) in favor over the two-norm (10). Clearly, 
the slew rate constraint (11) can be equivalently ex- 
pressed in terms of the three scalar conditions 

l -^i{t') oj max , i — 1,2,3. (64) 

Applying the findings of Lemma 1 in the Appendix 
to any one of the components in (64) we obtain the 
following result. 

Property 8 Consider the system (l)-(3) with the 
control law (12) and a slew rate constraint of the 
form (64) with fixed i £ {1,2,3} and co max > 0. 
Assume that the functional dependence of the gams 
k and c on their arguments u>, q, q$, and t is such 
that 

c>-k— (65) 

CO 4 

is satisfied whenever 

\coi\ A io max . ( 66 ) 

Then the slew rate constraint 

K-(t)| < ui max 

is guaranteed to remain satisfied for all times t >t 0 
if it is satisfied at the initial time to- 

Proof: Noting that 

LUi = -coJi - k qi (67) 

(see equation (13)), the statement of Property 8 fol- 
lows from applying Lemma 1 in the Appendix two 
times, namely to the constraints w; — co max < 0 and 
u)i oj max A 0. 

q.e.d. 

In complete analogy to Property 7 the statement 
of Property 8 can be strengthened by replacing the 
equality sign in (66) through an equality sign. Ex- 
plicitly, we have the following result. 
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Property 9 Consider ihe system (l)-(3) with the 
control law (12) and a slew rate constraint of the 
form (64) with fixed i £ {1,2,3} and co max > 0. 
Assume the gams k(u>, q, q$, t) and c(w, q, qq, t) 
are Lipschitz-bounded with respect to u>, q, 54 , and 
piecewise continuous with respect to t. Additionally, 
assume that k and c are such that 

c A — k —— ( 68 ) 

UJi 

is satisfied whenever 

\iOi \ = L0 max . ( 69 ) 

Then the slew rate constraint 

K-(t)| < ui max 

is guaranteed to remain satisfied for all times t > t o 
if it is satisfied at the initial time to- 

Proof: Follows from Lemma 2 in the Appendix. 

q.e.d. 

Properties 8 and 9 apply for each component of the 
three constraints i = 1, 2, 3 in (64) individually. Any 
combination of these constraints is satisfied if k and 
c are such that conditions (65), ( 66 ), in the case of 
Property 8 , and conditions ( 68 ), (69), in the case of 
Property 9, are satisfied for the appropriate combi- 
nation of indices i. 

In complete analogy to Remark 1 we have: 

Remark 3 Assume c and k are such that Prop- 
erty 8 (or Property 9) is satisfied. Let Ai > 0, 
A 2 > 0 be positive constants. Then, Property 8 (or 
Property 9) remains satisfied if we replace 

c —>■ c Ai 

k —>■ k X 2 


In complete analogy to the conclusions drawn from 
Remark 1 it follows from Remark 3 that, by picking 
c smaller (or k larger) than required by the condi- 
tions of Properties 8 or 9, the maximum attainable 
slew rate is reduced to values smaller than co max . 
In practice, this may represent an undesired per- 
formance degradation, and it may hence be a design 
goal to satisfy conditions (65) / ( 68 ) only marginally, 
i.e. with equality rather than strict inequality. From 
the same arguments as those used in the paragraphs 
following Remark 1 it follows, as in the previous sec- 
tion, that it is not sensible to blindly enforce condi- 
tions (65) / ( 68 ) with strict equality throughout. 

In analogy to Remark 2, a convenient simplifica- 
tion of the conditions of Properties 6 and 7 can be 
given as follows. 


Remark 4 The conditions (65), (66) of Property 8 
(conditions (68), (69) of Property 9) are guaranteed 
to be satisfied if the functional dependence of k and 
c on their arguments u>, q, q 4 , and t is such that 

c= 1*1 J«iL. (70) 

whenever |w;| > uj max (whenever |w;| = ui max ). 

Summary and Conclusions 

We have analyzed a nonlinear feedback approach 
to the problem of reorienting a rigid spacecraft from 
arbitrary initial conditions to prescribed final con- 
ditions with zero angular velocity. The control 
law, which is derived from feedback linearization, 
is based on quaternion feedback and leaves the user 
to choose two gains as functions of position, angu- 
lar rate, and time. Under mild conditions on these 
gains it is shown that the control law has the follow- 
ing properties: 1 ) rest-to-rest boundary conditions 
result in a principal axis rotation, 2 ) for arbitrary 
initial conditions, the rotational motion perpendic- 
ular to the principal axis decays exponentially, and 
3) the prescribed final states are guaranteed to be 
reached asymptotically in time. Additionally, the 
paper analyzes slew rate constraints, both, in terms 
of the two-norm and the infinity-norm. For both 
cases, conditions on the controller gains are derived 
that guarantee satisfaction of the slew rate limit at 
all times. 

In summary, the paper provides the design engi- 
neer with a general nonlinear control law with ap- 
pealing closed-loop properties. Global asymptotic 
stability and satisfaction of slew rate constraints can 
be guaranteed by observing the simple rules on the 
controller gain design derived in this paper. 
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Appendix 

Lemma 1 Consider the initial value problem 

x = f(x,t), x(t 0 ) = x 0 (71) 

and the state constraint 

g(x,t)< 0 (72) 
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Assume that the functional dependence of the func- 
tions f and g on their arguments x and t is such 


that 


dg(x,t) 


f(x,t) 


dg(x,t) 


< 0 


dx ’ 7 dt 

is satisfied for all x and t that are such that 


(73) 


g(x,t)>0. (74) 


Then the constraint (72) is guaranteed to remain 
satisfied for all times t > t o along the solution 
of (71) if (72) is satisfied at the initial time to, i.e., 
if g{x 0 ,to) < 0. 


Proof: We have 


dt 


g(x,t) 


dg(x,t) dg(x, t) 

8x 8t 


(75) 


Hence, conditions (73), (74) imply that ^ < 0 
whenever g > 0. Hence, with g(xo,to) < 0, it is 
clear that (72) remains satisfied for all times t > to- 


q.e.d. 


Intuitively, it can be expected that Lemma 1 re- 
mains satisfied if the inequality sign in (74) is re- 
placed by an equality sign. It can be shown that this 
is in fact correct if the functions f(x,t) and g(x,t) 
are sufficiently “well-behaved” . Explicitly, we have 
the following result. 

Lemma 2 Consider the initial value problem (71) 
and the state constraint (72). Assume that g is dif- 
ferentiable and that f(x,t), , and are 

Lipschitz-bounded with respect to x, and piecewise 
continuous with respect to t. Additionally, assume 
that the functional dependence of the functions f 
and g on their arguments x and t is such that 


dg(x,t) 

dx 


f(x,t) 


dg(x,t) < Q 
dt ~ 


is satisfied for all x and t that are such that 


(76) 


g(x,t) = 0. (77) 

Then the constraint (72) is guaranteed to remain 
satisfied for all times t > t o along the solution 
of (71) if (72) is satisfied at the initial time to, i.e., 
if g(xo,t 0 ) < 0. 

Proof: Define 


= ^/(*C) + ^ (78) 

and consider the initial value problem 


X 


l(x,t) 


x(to) 


Xo 

. y . 


h(x, t) 

5 

_ y{t o) . 


_ g(xo,t 0 ) 


(79) 


where 


h(x, y,t) 


h(x, t) if y < 0 
0 if y > 0 


(80) 


Let the solution of (79) be denoted by superscript 
hat. Using Lemma 1 it can be shown that the 
solution of the initial value problem (79) satisfies 
the constraint y < 0, i.e. we have y(t) < 0, if 
g(xo,to) < 0. Hence x(t), y(t) also furnish a solu- 
tion to the initial value problem 


X 


f(x,t) 


x(to) 


Xo 

. y . 


h(x, t) 

5 

_ y(to) 


_ g(xo,t 0 ) 


9 (81) 

if g(xo,to) < 0. Clearly, the functions /, and 
^ in (81) are such that it can be guaranteed that 
the solution to (81) is unique (see the Appendix I 
of [13]). Together with the results obtained above, 
this implies that the unique solution to (81) satisfies 
y(t) < 0. Noting that the j/-component in (81) can 
be integrated analytically to yield y(t) = g(x(t),t) 
proves the statement of the lemma. 

q.e.d. 


Remark 5 The statement of Lemma 2 can be 
strengthened by replacing “piecewise continuous with 
respect to time t” through “measurable with respect 
to time t. ” 


Remark 6 The assumption of Lipschitz bounded- 
ness in Lemma 2 is, in fact, necessary and cannot 
be weakened. 


For example, consider the initial value problem 
x = \/ , x(to) = xo, with to = xo = —1 (82) 

and the constraint 

x < 0 (83) 

The solution x(t) of (82) obviously satisfies x(t o) < 0 
and the right-hand side of the differential equations 
in (82) are such that x < 0 whenever x = 0. How- 
ever, it can be easily verified that 

x(t) = t 3 

furnishes a solution to (82), which obviously violates 
x < 0 for t > 0. In fact, the right-hand side of the 
differential equation in (82) is not Lipschitz bounded 
at x = 0, and it can be verified that 

ft 3 for t < 0 
~ \ 0 for t > 0 

represents an alternative solution to (82). 
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